Couple of points here:

• Static cases like axial stretching and timoshenko are okay. The analytical results are recovered
• Catenary case blows up with the current parameter. Might need to adjust the damping constant
• Snake and flagella cases are broken due to some issue in base_system (e.g. run the snake case to replicate the issue). It does not seem to be related to this PR.
• The new implementation is okay, but it erases the original behavior and may break user code after the update is released. I thought the plan was to keep the original behavior but add a deprecation warning?

Catenary case blows up with the current parameter. Might need to adjust the damping constant

@Ali-7800 Do you remember how you tuned the damping coefficient for Catenary case before?

• Snake and flagella cases are broken due to some issue in base_system (e.g. run the snake case to replicate the issue). It does not seem to be related to this PR.

Those cases will be fixed with #367

• The new implementation is okay, but it erases the original behavior and may break user code after the update is released. I thought the plan was to keep the original behavior but add a deprecation warning?

We can do that if we consider this change as an update. I think we can make a hotfix if we consider this change a bug fix. @armantekinalp opinion?